Gauging a quantum heat bath with dissipative Landau-Zener transitions 
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We calculate the exact Landau-Zener transitions probabilities for a qubit with arbitrary linear 
coupling to a bath at zero temperature. The final quantum state exhibits a peculiar entanglement 
between the qubit and the bath. In the special case of a diagonal coupling, the bath does not 
influence the transition probability, whatever the speed of the Landau-Zener sweep. It is proposed 
to use Landau-Zener transitions to determine both the reorganization energy and the integrated 
spectral density of the bath. Possible applications include circuit QED and molecular nanomagnets. 
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Quite a number of quantum two-state systems are 
presently tested as candidate qubits, the units of quan- 
tum information. Good qubits are well isolated from 
their environment, but easy to manipulate. This some- 
what conflicting requirement has spurred renewed inter- 
est in the dynamics of qubits coupled to an environment 
or heat bath 0,0- Qubits can be seen as bath detectors. 
For example, static qubits probe via their decay rates 
the bath spectral density at their transition frequencies. 
In solid-state environments these rates can be strongly 
frequency- and sample-dependent. We discuss an in this 
respect superior 'bath detection mode' of the qubit. 

One way of changing the state of a two-level sys- 
tem involves the forced crossing of its diabatic energies. 
For constant level-crossing speed this is known as the 
Landau-Zener (LZ) problem 0, 0,0, which for a two- 
level system can be solved exactly jj, y, M ■ This is no 
longer the case when taking its environment into account 
0, Il0j | that may cause thermal excitation and quan- 
tum tunneling. In the low-temperature tunneling regime, 
analytical estimates for transition probabilities exist only 
for very fast and very slow sweeps, and the literature is 
not unanimous about the latter limit 0, S 0, Hoi Hl| . 

In another line of research, incited by the paper |l2| . 
it was recently proven that LZ transition probabilities 
can be calculated exactly for some many- level systems as 
well, although only for some initial states [3E1E1E3- 
In this Letter, we extend the analysis to quantum dissipa- 
tive systems and study LZ transitions in spin-boson prob- 
lems in a new way. First we calculate zero-temperature 
LZ transition probabilities exactly. Second, by making 
use of this exact dependence on the bath parameters, we 
propose to gauge the dissipative environment of a qubit 
by performing LZ sweeps. One advantage of this bath 
detection mode of the qubit is that effects of spiky bath 
spectral densities are averaged out in every sweep. 

Driven spin-boson model. — Consider LZ transitions in 
a qubit coupled to a bath of N quantum harmonic oscil- 



lators, as described by the Hamiltonian 



JV 



H(t) = ycr^ + ^a x + Kfybtbj + H QO , (1) 



4=1 



with the qubit-oscillator coupling 



4=1 



JV 

y [ cos °j (T z + sin 6j(T x ] {bj + b]). (2) 



The energy difference between the diabatic qubit states 
changes linearly in time as vt (with level-crossing speed 
v > 0) and their intrinsic interaction amplitude is A. The 
cr x . z are Pauli operators. The first two terms of define 
the standard Landau-Zener problem for an isolated two- 
level system. The N harmonic oscillators in JQ) can have 
different frequencies £lj > 0, qubit-oscillator couplings 
7j, and interaction angles 6j. The oscillators affect the 
qubit (i) by changing its energies via the diagonal cou- 
pling oc jj cos 6 j cr z , and (ii) by inducing transitions be- 
tween its levels via the transverse coupling oc "fj sin 6j(T x . 

In order to calculate the probability that the qubit 
state flips due to the LZ sweep, we will work in an in- 
teraction picture and split the Hamiltonian QJ into an 

interaction V = cr x [A/2 + J2jLi (74 / 2 ) sin #4 ( b j + b })] in- 
volving bit flips as described by cr x = |f)(j| + |I)(T|, 
and the (bit-flip-)frce Hamiltonian Ho(t). A polaron 
transformation 2] diagonalizes Ho(t) in terms of qubit- 
state dependent shifted oscillators, giving Ho(t) the form 

Ea= T >X«l{±f + Y.U m A± h i± ~ E \ with cre - 

ation and annihilation operators for the shifted oscil- 



lators 6^2 = b^' ± Qj/2 with Q 3 = {■yj/W.j) cos^-; 
the + corresponds to the qubit state |f) and the — 
to j!). With bath oscillators shifting in this way, 
the reorganization energy Q gained by the system 

has the same value E = J2j=i(lj/^^j) cos2 ^j f° r 
both qubit states. Eigenstates of the shifted oscilla- 
tors are labelled as |n±), where the N components 



^(t) 



2 



rij of the vector n are single-oscillator excitation num- 
bers. The free time-evolution operator Uo(t) can be writ- 
ten as E s=T ,i \s)(s\e Tivt2/ih J2n l n ±)( n ±l e - in - nt+iEt ^ , 
where the inner product of the vectors n and 11 = 
(Oi, Q2, ■ ■ ■ ,&n) shows up. We next define the 
interaction-picture Hamiltonian as H (t) = U^t) V Uo(t). 
In order to bring H(t) into a useful form, we write the 
oscillator operators (bj + &t) in V as (bj± + b) ± T Qj)- 
We choose to associate the "-(-"-oscillators with the term 
||) (|| of cr x , and the "—"-operators with the other term. 
We then write bj± = J2 n \f™j\( n ~ lj)±)( n ±|j where 
lj is the unit vector with jth component equal to 1, 



and likewise bj ± = J2 n \J n i + l|(n + lj)±)(n±|. In this 
form, the 'bra's' of shifted oscillator states in V corre- 
spond to the 'bra's' of the qubit, while for the 'kets' this 
must still be arranged by using the completeness rela- 
tion |m±) = ^ n |n =F )(n =F |m±). The interaction-picture 
Hamiltonian then becomes 

H(t) = ^ e i(m - n) - m {^mnlT>a| e iu * 2/2 ^|m + >(n_| 

III . II 

+^ mn |l)<T| e ^ t2/2?l ®|m_)<n + |}, (3) 
involving the two infinite-dimensional matrices W with 



JV 



W ± = 

rim 



A±^ 7j g j sin( 



N 

!m±|n T ) +^7 J sin^{^/nJ(m ± |(n- lj) T ) + y/nj + l(m±|(n- 



>}• (4) 



Transition probability. — We focus on the situation that at time t — — oo the system starts in its ground state 
\tl>(— oo)) = |t) |0+). We are now interested in the survival probability P|-_,i-(oo) of the initial state |t) of the qubit. 
-Pj^l(oo) equals the square of the norm of the projected final oscillator state ('\\ip(oo)}, which can be written as 
(1\U(oo, — oo)||) |0 + ), where U(t2,t\) — T exp[— i\/fi) dTH(r)} is the time evolution operator in the interaction 
picture. In a time-ordered expansion of [/(oo, — oo), only the even powers of H(t) will contribute to Pf^f (oo). As is 
well known, the (2fc)th-order term in the expansion involves a (2fc)-fold time integral with variables t2k > fak-i > ■ ■ ■ > 
t2 > t\ in the interval (— oo, oo). It is advantageous to make the variable transformation 6] x q — t\ +X}l=i {fat+i — fat) 
for 2 < q < k, x\ = ti, and y q = t q — biq—i. We label the iV-oscillator state after £ interactions as . For brevity, we 
define the frequencies we = (n^ — n^ -1 )) • ft and wi = • fl. The perturbation series for (||^(oo)) then becomes 



^ ^ POO POO POO POO POO 

YsTmxIk- da; i / d ^2.../ dx k d Vl ... dy k ^»)„(»-i)^-i)„(«- 2 ) ••■ w Zm a M W *n 

fe =0 ^ I J -°° Jx i Jx >*-i J ° J ° „(2*0 ... n (D 



x exp 



k ( h l E ~ x \ k ( f_1 N 

^ xtyi + - (w 2 £ + w 2 e-i)x e + -y\ + ye'Y / yq\ w 2 eye + {w 2 e +W2i-i)^2y q 



i=i 



9=1 



9=1 



I (2fc)\ 

\n\ '). 



(5) 



It would be convenient to symmetrize the a^-integrals 
at this point, but in general the integrand is not sym- 
metric under permutation of the X£. When transforming 
the Xi to new variables si = x\ and se = xg — xe-i for 
I = 2, 3, . . . , k, one finds that the /_ dsi -integral yields 

the delta-function (2^/v)8(J2e=i Vt + Ji{ w 2t + w 2e-i) )■ 
Now, since the initial state is ||) |0+), the sum over the 
we in this delta function equals n^ 2 ^ ■ ft, which evi- 
dently is > 0. Likewise, the variables ye are all pos- 
itive or zero. Therefore, the delta-function can only 
"click" in the subspace |/i = = ... = y2k = 
and will do so only if the vector n( 2fc ) vanishes. The 
physical meaning of the latter statement is discussed be- 
low. Performing the other s^-integrals is cumbersome, 
so we quickly return to the integrals in Eq. but this 
time armed with the knowledge that only the subspace 
2/i = ... = y2k = contributes. Since within this sub- 



space the integrand is symmetric in the variables xe, it 
is correct to symmetrize the a^-integrals in Eq. JSJ), i.e. 
we can replace them by (1/fc!) da^i . . . dxk- Af- 
ter performing these standard integrals, the ^-integrals 
can be evaluated as well. For example, the yi-integral 
J dyi S(yi+hw 2 /v+hwi/v) vanishes unless (u>2+wi) = 
rr 2 ' -ft vanishes, in which case it is indeed the "subspace" 
defined by y\ = that makes this integral equal to \. 

From the time integrals we find the following selec- 
tion rule: when starting in the ground state ||) |0 + ), the 
only (2fc)th-order processes contributing to the survival 
probability Pf- >-T (°°) are those with n^ 2 ) = n^ 4 ) = . . . = 
n (2fe) _ o Hence the oscillators will end up in their ini- 
tial state |0+) in case the qubit ends up in ||). This 
stri king result agrees with the so-called no-go theorem 
[l2l ll3L H3. Il5j , which we extended here to spin-boson 
problems. The time integrals do not forbid occupation 



3 



of the states |T)l n + 7^ 0+) at intermediate times, nor do 
they restrict the intermediate oscillator states |ni _1 ^), 
but further restrictions may originate from vanishing ma- 
trix elements W^ n , see Eq. (QJ. Only when the qubit 
ends up in ||) can the LZ driving dissipate energy into 
the bath. Hence qubit and bath end up entangled. 

In line with the selection rule, we find that (tlV'(oo)) 
simplifies into exp(— ttW 2 /4Tiv) |0 + ), where the parame- 
ter W 2 =Yln Won^no 1S still to be determined by using 
Eq. We finally obtain our central result: The exact 
Landau-Zener transition probability for a qubit arbitrar- 
ily coupled to an oscillator bath at T = is 

P T ^(oc) = f - P t .,(oo) = 1 - e -^ 2 /«, (6) 

where the parameter W 2 is given by 

n / sin#, ;cos#,'7?\ 2 J-^ „ „ 

^ 2 = ( A -E nn 3 ) + I> 2 ^- (7) 

V 3=1 3 7 3 = 1 

By introducing the spectral density J(w) = 
^ 3=1 (27j/?i) 2 (5(cj — flj) and assuming that oscilla- 
tors with equal frequencies to have equal coupling angles 
6(u>), Eq. becomes 

W 2 = J sm[9(u)}cos[9{uj)}J{uj)/ujj 

h 2 r°° 

+_ / du sin 2 [6 (w)] J (cA (8) 

No specific form of J(w) or 9(uj) is presupposed, but for 
the limit N — > 00 describing a continuum of oscillators, 
the two integrals in Eq. (jSJ have to be finite. 

Transverse coupling (9j — tt/2). — We first focus on 
purely transverse bath coupling and obtain W 2 = A 2 +S, 
where S = J2f=i l] = i^/i^)) Jo°° dwJ(w) is the inte- 
grated spectral density. The jth oscillator reduces the 
final survival probability by a factor exp[— it^ 2 /{2hv)], 
independent of the other oscillators or of the value of flj . 
For A = 0, the bath ends in |0+) when the qubit ends in 
|t), whereas the bath contains an odd number of bosons 
when the qubit ends in ||). Qubit and bath therefore 
end up fully entangled, in the sense that the qubit coher- 
ence vanishes after tracing out the bath. Nonzero final 
qubit coherence requires A ^ 0. Notice also the interest- 
ing phenomenon of bath- assisted adiabatic following: For 
large bath coupling, i.e. S 2h~v/ir, the qubit ends up 
in the initially unpopulated state, even in the absence of 
an intrinsic interaction A. 

Diagonal coupling (9j = 0). — In solid-state environ- 
ments, decoherence often occurs much faster than the 
bath-induced relaxation so that the latter is neglected. 
Hence the "standard model" for LZ transitions in dissi- 
pative environments 0, H, is the single qubit di- 
agonally coupled to an oscillator bath, which is obtained 
from our Hamiltonian by setting all 9j =0. By doing 



the same in Eq. (JTJ, we find that for a qubit diagonally 
coupled to a bath at T = 0, the Landau-Zener transition 
probability Pf^j(oo) = 1 — cxp[— irA 2 /(2hv)], which is 
the well-known result for an isolated qubit! Thus, al- 
though the bath coupling does not commute with the 
qubit Hamiltonian, there is no bath dependence in the 
transition probability, neither via qubit-bath couplings 
7j nor via the oscillator frequencies Qj, no matter how 
fast the LZ sweep is performed. We discuss this below. 

Gauging quantum dissipation. — Next assume that all 
oscillators couple to the qubit with the same angle 9j = 
8. This may follow from microscopic considerations or 
may be engineered. For example, when all oscillators 
couple diagonally to the qubit in a basis described by 
Pauli matrices T x y z , then an experiment where the qubit 
is driven via (vt/2)a z = (vt/2)(cos0T z — sin 6t x ) while 
(A/2)cr x = (A/2) (sin 8t z + cos 9t x ) is kept constant, is 
described by our Hamiltonian with all angles 8j = 9. 
The reorganization energy then becomes E = Eq cos 2 9, 
where E = ^=1^/(^3) = W ^) J™ dto J {lu) / lo 
occurs for diagonal coupling. From Eq. (JJJ) we deduce 

W 2 = (A - E sin 9 cos 9) 2 + S sin 2 9. (9) 

Unlike for transverse coupling, the oscillators do not af- 
fect -Pf^T (00) independently. Relaxation oc sin^, dephas- 
ing cx cos 9, and intrinsic interaction act together. In the 
continuum limit N — > 00, Eq. holds for spectral den- 
sities J(uj) that give rise to finite reorganization energies 
E and finite integrated spectral densities S. 

As a main application of our calculations, we propose 
to gauge the dissipative environment of a qubit via LZ 
transitions. By gauging we mean in short: the measure- 
ment of Eq and S and the subsequent parameter fixing 
of appropriate model spectral densities. More in detail, 
we propose to perform LZ sweeps and to determine tran- 
sition probabilities for several fixed values of the tunable 
intrinsic interaction A. Notice that H^ 2 (A) as a function 
of A is a parabola in Eq. © . Figure depicts corre- 
sponding final transition probabilities for several coupling 
angles. Diagonal (9j = 0) and transverse (9j = tt/2) 
coupling are limiting cases which have in common that 
IF 2 (A) and hence Pf_>|(oo) assume their minima for 
A = 0. Neither limiting case is suitable for gauging 
Eq. Since for intermediate cases < 9 < ir/2 we find 
a minimal transition probability for a nonzero internal 
interaction A m j n = EqSih9 cos9 7 one can determine the 
reorganization energy Eq by measuring A m ; n . Moreover, 
the integrated spectral density S can then be identified 
as S = VF 2 (A m i n )/ sin 2 9. As a consistency test, one can 
check whether W 2 (A = 0) = sin 2 9(S + E 2 cos 2 9) holds. 
Consistency can be further tested by changing the basis 
in which the LZ sweep is performed: the values of Eq 
and S* should come out independent of the angle 9. If 
not, then 9(uj) is not constant and Eq. (JHJ) applies. 

With the values of Eq and S thus determined, one can 
fix parameters in appropriate model spectral densities. 



4 



1 

0.75 
0.5 
0.25 










S / / ' 
/ f / 

/ / / — " o 




// — ■ 6> = tt/8 


~ x V 
\ / . 
\ / 


- - 6» = tt/4 




(9 = 37r/8 


A 

y , i 


— <9 = 7r/2 




0.5 



1 1.5 2 
A/Vhv 



2.5 



FIG. 1: (Color online) Final transition probability Pf^x(oo) 
as a function of intrinsic interaction A, for several values of the 
coupling angle 8. Parameters: Eo = 2.0Vhv and S = 0.5hv . 



For example, suppose J(oj) has the known form J(uj) — 
auj s e^^^ c with a power s > given by the physical 
nature of the environment , but with unkown strength 
a and cutoff frequency lo c . (s = 1 is the Ohmic case and 
for s < 0, E would diverge.) Then E = ahuj^T (s) / (Aw) , 
in terms of the gamma function, and S = ah uj^ +1 T(s + 
l)/(47r). Hence LZ gauging fixes huj c = S/(E s) and a. 
Other spectral densities may depend on more parameters. 
LZ gauging fixes two of them. 

The efficient use of our gauging scenario requires the 
following: For i-j^(oo) to change considerably, one 
chooses v such that (nS/2hv) sin 2 ~ 1 and varies A 
on the scale of \/Ssm9. Then Eq can be measured accu- 
rately if Eo cos 9 is not much smaller than \/~S, or \J a/ Air 
not much smaller than unity for an Ohmic bath. So 
Pf^j(oo) is robust under dephasing, even if dephasing 
is much faster than relaxation (9 <C 1), unless the qubit- 
bath coupling is strong. This robustness has found its 
use in experiments on molecular nanomagnets [ItI llif . 

Discussion and conclusions. — We presented in 
Eqs. © and J7|) exact Landau-Zener transition probabil- 
ities for a qubit with arbitrary linear coupling to a zero- 
temperature bath. We found that qubit and bath end 
up entangled. Our results apply to experiments where 
the initial and final qubit energies are off-resonant with 
relevant bath frequencies, for example in circuit QED 
|l9l lifH , where qubit energies can be varied over a broad 
range and spectral densities are peaked 0, 122^ . In- 
deed, our predictions for transverse coupling generalize 
our detailed study of LZ sweeps in circuit QED [16| to 
more realistic situations where peaked spectral densities 
have nonzero widths. Other applications include tunable 
atoms in optical cavities and in photonic crystals. 

For diagonal coupling we find that the transition prob- 
ability does not depend on the bath at all. Our exact re- 
sult settles a long-standing discussion 0, IE 0, H EI , at 
least for zero temperature. It corroborates and interpo- 



lates between what Ao and Rammer gj, |9| found for the 
fast-passage limit A 2 /hv <C 1 and for the opposite adia- 
batic limit A 2 ITiv <C 1, which both have been confirmed 
numerically [l(j|. For fast LZ sweeps, the bath clearly has 
no time to affect the transition, but the absence of any 
bath influence also for slower sweeps is a highly nontrivial 
property of this standard model. 

A qubit undergoing a LZ sweep measures 'global' prop- 
erties of the bath, namely the frequency-integrated spec- 
tral density, and for strong coupling the reorganization 
energy. Sample-dependent spikes in spectral densities are 
therefore averaged out, even in a single sweep, so that 
model parameters can be determined. We therefore pro- 
pose to employ LZ transitions for a valuable gauging of 
the dissipative environment of a tunable qubit. 
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